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Measurement of halo properties with weak lensing shear and flexion 
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ABSTRACT 

We constrain properties of cluster haloes by performing likelihood analysis using lensing 
shear and flexion data. We test our analysis using two mock cluster haloes: an isothermal 
ellipsoid (SIE) model and a more realistic elliptical Navarro-Frenk- White (eNFW) model. 
For both haloes, we find that flexion is more sensitive to the halo ellipticity than shear. The 
introduction of flexion information significantly improves the constraints on halo ellipticity, 
orientation and mass. We also point out that there is a degeneracy between the mass and the 
ellipticity of SIE models in the lensing signal. 
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1 INTRODUCTION 

The properties of galaxy and cluster haloes are of great interest in 
cosmology, and can be powerful tests of the cosmological paradigm 
and the nature of dark matter. Two important parameters that de- 
scribe a dark matter halo are its mass and shape, which are related 
to many physical processes, such as the growth and merging history 
(Kauffmann et al. 1993; Springel et al. 2005). 

Models of dark matter haloes beyond the spherical approxi- 
mation are favored by many numerical simulations (Jing & Suto 
2002; Springel et al. 2004; Kazantzidis et al. 2004; Allgood et al. 
2006) and observations (Reblinsky 2000; Lee & Suto 2004; De 
Filippis et al. 2005; Sereno et al. 2006; Wang et al. 2010). Fur- 
thermore, numerical simulations with different assumptions predict 
different properties of dark matter haloes (e.g. Bullock 2002; Bailin 
& Steinmetz 2005; Wang & White 2007). Current models based on 
N-body simulations, semi-analytic models or hydrodynamic simu- 
lations can predict several halo properties, but several ingredients 
of these models remain uncertain. A precise understanding of halo 
properties such as mass and ellipticity is important to confirm and 
improve the existing models of galaxy formation and probe the 
physical nature of dark matter. 

Gravitational lensing is a powerful tool to study mass distri- 
butions, independent of the nature or dynamical state of the mat- 
ter (see Mellier 1999; Schneider 2006; Munshi et al. 2008, for re- 
views). Galaxy-Galaxy Lensing (GGL) is concerned with the mass 
associated with galaxies and dark matter haloes in which galaxies 
reside (Tyson et al. 1984; Brainerd et al. 1996; Hudson et al. 1998; 
Hoekstra et al. 2004; Sheldon et al. 2004; Mandelbaum et al. 2006). 
The distortion caused by a single galaxy cannot be detected, but the 
statistics of many foreground-background pairs yield a detectable 
signal for a population of galaxies. Brainerd et al. (1996) discov- 
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ered a significant GGL shear signal. Schneider & Rix (1997) de- 
veloped a maximum likelihood analysis that can constrain the halo 
properties of the lens galaxy populations through GGL, allowing to 
estimate the mean velocity dispersion and the characteristic scale 
for a non-singular isothermal sphere halo model. 

Flexion as the gradient of the projected mass density, is sen- 
sitive to the small-scale variations of mass distributions (Gold- 
berg & Natarajan 2002; Goldberg & Bacon 2005; Bacon et al. 
2006). Different techniques have been developed to measure flex- 
ion (see Irwin & Shmakova 2006; Okura et al. 2007; Schneider & 
Er 2008; Fluke & Lasky 2011, for examples). Recently, Velander 
et al. (2011) applied the shapelets technique on the COSMOS sur- 
vey, and Cain et al. (201 1) introduced a new method (so-called an- 
alytic image model) to analyse lensing flexion images. It has been 
noted that flexion can contribute to cosmology in several aspects, 
such as exploring the mass distribution of dark matter haloes of 
galaxies and clusters, especially substructures (Leonard et al. 2009; 
Bacon et al. 2010; Er et al. 2010). Hilbert et al. (201 1) also propose 
to reduce the distance measurement errors of standard candles us- 
ing lensing shear and flexion maps. 

In Hawken & Bridle (2009); Er & Schneider (201 1); Er et al. 
(2011), the ellipticity of a galaxy halo has been studied with flex- 
ion. It was found that the constraints from flexion are tighter than 
those from shear. In Goldberg & Bacon (2005), Galaxy-Galaxy 
lensing Flexion (GGF) has been studied using the distribution func- 
tion of orientations between the line connecting foreground and 
background pairs and the flexion of the background galaxies. A 
GGF signal has been detected by Leonard et al. (2007) using im- 
ages taken by HST ACS in the cluster Abell 1689. Moreover, com- 
bining shear and flexion information provides tighter constraints on 
halo properties by studying mass distribution on different scales (Er 
et al. 2010; Shapiro et al. 2010; Hilbert et al. 201 1). 

In this paper, we combine shear and flexion data to constrain 
the properties of dark matter haloes. The tangential shear is mainly 
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sensitive to the mass of the halo, whereas flexion is sensitive to the 
halo ellipticity. Moreover the usable number density of flexion data 
is relatively low, since the flexion signal drops faster than the shear 
signal with the angular distance to the centre of the halo. It is thus 
not sufficient to constrain the halo properties with shear or flexion 
alone. Therefore, we propose to take advantage of both shear and 
flexion in our analysis. In particular, we use the angular positions, 
tangential shear and flexion of galaxies in our likelihood functions. 
A singular isothermal ellipsoid model and an elliptical NFW model 
for a galaxy or cluster halo with 3 parameters (mass, ellipticity and 
orientation) are adopted in this paper. In Sec. 2, we recall the basic 
lensing equations. Our likelihood function is introduced in Sec. 3. 
We perform numerical tests of our method and results in Sec. 4 
and present our conclusions in Sec. 5. Throughout this paper, we 
adopt a ACDM model with fl A = 0.75, f2 m = 0.25, and a Hubble 
constant Ho = 73 kms -1 Mpc~\ 



2 LENSING BASICS 

The formalism described here can be found in Schneider & Er 
(2008); Er & Schneider (201 1). The weak lensing shear and flexion 
are conveniently described using a complex formalism. We adopt 
the thin lens approximation, assuming that the lensing mass dis- 
tribution is projected onto a single lens plane. The dimensionless 
projected mass density can be written as k(6) = E(0) /E cr , where 
is the vector of (angular) position coordinates, E(0) is the pro- 
jected mass density and E cr is the critical density, given by 

1 D s (oo) 



4ttG D d D d . 



(1) 



for a fiducial source located at a redshift z — > oo. Here D s (oo), D d 
and Dd.oc are the angular diameter distances between the observer 
and the source, the observer and the lens and between the lens and 
the source, respectively. 

The first order image distortion induced by gravitational lens- 
ing is the shear 7, which transforms a circular source into an ellip- 
tical one. The second order effect, called flexion, is described by 
two parameters: the spin-1 flexion, which is the complex derivative 
of K 



^ = VcK =^ +i ^' (2) 
and the spin-3 flexion, which is the complex derivative of 7 



Q = V c7 . 



(3) 



For a source at redshift z B and a lens at redshift z d , a 'cosmo- 
logical weight' function must be introduced: 



Z(z s 
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(4) 



where Dds and D B are the angular diameter distances between the 
lens and the source, and the observer and the source. H(z s — Zd) 
is the Heaviside step function to ensure that the source redshift is 
higher than the lens redshift. The first and second-order lensing 
effects scale with the source redshift as 



k(z s ) — Z(z s )k, 
F{z B ) = Z(z s )T, 



70s) = Z(z s )j, 

g(z s ) = z{z s )g. 



(5) 



In GGL we express the shear with respect to a foreground 
halo. This defines the tangential shear 



where 71, 72 are respectively the real and the imaginary compo- 
nents of shear, and ip is the polar angle with respect to the vector 
connecting the background and foreground galaxies. T\ {Qi) and 
Ti (G2) are the real and imaginary components of the spin-1 (spin- 
3) flexion. 

In this paper, 7t, J-i, Ti, Qi and Q2 represent the values cal- 
culated from our model, while e, J"f bs , T% hs , Gf s and G% bs are 
the observables, which will be introduced in the next section. 



3 METHODOLOGY 

We apply the Bayesian framework to study galaxy-galaxy lensing 
simulated maps, in order to estimate physical parameters of a fore- 
ground dark matter halo. The maps contain tangential shear and 
two flexion components at the positions of background galaxies. 
The values of these observables reflect the intrinsic shapes of the 
background images, slightly modified by lensing due to the gravi- 
tational potential of a foreground halo. 

We assume the shear and flexion components can be measured 
with unbiased estimators that are linear combinations of the lens- 
ing and intrinsic contributions to an image shape. The ellipticity of 
background galaxies e is such an estimator. Analogously, we as- 
sume that the estimators of flexion, J? bs , J^* 8 , G? hs and G% hs are 
derived from higher-order brightness moments of the background 
images. We further assume that the noise of such estimators is un- 
correlated due to the distributions of intrinsic shapes, and neglect 
other contributions, such as sample variance. 

The intrinsic ellipticity distribution can be described by a 
Gaussian probability density distribution with zero mean and stan- 
dard deviation a e « 0.3 (Brainerd et al. 1996). For flexion, we 
use the preliminary studies of intrinsic flexion by Goldberg & 
Bacon (2005) and Goldberg & Leonard (2007), where scatter of 
a? — 0.03 per arcsecond and ag = 0.04 per arcsecond were found 
for spin-1 and spin-3 flexion, respectively. In this paper, we also as- 
sume that the distribution of each intrinsic flexion noise component 
riFj — T 

zero mean and a T = 0.03/", ag = 0.04/" 



- Tj, nGj — Qf° s — Gj = 1) 2) is a Gaussian with 



P(n Fj ) 



1 



2%ag 



exp 



exp 



2ol 



(7) 



(8) 



Constraints on the foreground haloes parameters p, which will 
be introduced in the next section, are inferred by evaluating their 
likelihood functions. The likelihood of a model is the conditional 
probability of the data given the model. Following the previous dis- 
cussion, we define three likelihood functions: 



L F {Ji{p)) = 
La{Gi{p)) = 



2-7T a, 
1 



■exp 



-(e» - 7O 2 



2-7T ffjr 

1 



exp 



2ai 



2a% 



2% ag 



exp 



(g° bs -Qif 



2a% 



(9) 



(10) 



(11) 



7t 



-71 cos 2ip — 72 sin 2ip, 



(6) 



The subscript i refers to the i* background galaxy. Notice that not 
all images have a measurable shape and usually flexion is measur- 
able only for a subset of all the images. For this reason, a selection 
must be applied to the data. The details of the mock data produced 
are given in the next section. 

The galaxy shapes are assumed to be independent, i.e., in this 
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Figure 1. The coordinate systems and relative angles used. The star repre- 
sents the location of one background galaxy. 



with p(0i , 02) defined by 



(14) 



The halo produces the following shear and flexion fields for a 
source at z s — oo: 



0e 

•T(gi,fl£) 

Q(o'i,e' 2 ) 
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(15) 



(16) 
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+ 



9 2 p 3 



^ 60^02 - 20 2 3 0f0 2 -0| + 2/ 2 0^0 2 ; 
4 p + 2 p 3 ^ 



where = \/0f + 0| and / is the axial ratio. The shear and flexion 
will be scaled to different redshifts for other values of z s , according 
to Eq. (5). 



study we assume uncorrelated spatial noise, and no systematic in- 
trinsic shape correlations or spurious correlations from PSF resid- 
uals. Each of the three likelihoods can thus be multiplied over all 
galaxy pairs with measurable background galaxy shape informa- 
tion. Furthermore, if we assume that the measurements of tangen- 
tial shear and flexion are independent at each galaxy position, the 
likelihood can be written as 



Y[{Le{0i)) 



n(M*')) 



\[{L G {ei)) 



(12) 



Recently Viola et al. (201 1) pointed out a correlation between shear 
and flexion noises. Although the correlation can be reduced, a pre- 
cise covariance treatment of shear and flexion will require more 
detailed empirical knowledge of the flexion noise and is beyond 
the scope of this paper. 



4 ANALYSIS AND RESULTS 

In this section, we first use the analytic singular isothermal ellip- 
soid (SIE) model to illustrate the results before we study the more 
realistic elliptical Navarro-Frenk- White (eNFW) profile. 



4.1 SIE model 

We adopt an SIE model for the foreground lens. The SIE model 
includes 3 parameters: the Einstein radius 0e, which defines the 
scale of the lens and is related to the mass or velocity dispersion 
(seeEq. 18); the halo ellipticity, defined by e = (9 a -6 b )/(6 a + 6 b ) 
(or equivalently by the axial ratio / = 9 b /9 a ), where 9 a , 8 b are the 
major and minor axes; and the halo orientation (f>o. 

The lensing properties of the halo, such as shear and flexion, 
are calculated at the background galaxies positions (0i,02). The 
image reference frame relates to the halo reference frame through 
the rotation 0i = 6[coscj>o + 02sin0o, 02 = — 0isin<^o + 
2 cos 4>o , where cj>o is the halo orientation (see Fig. 1 for an il- 
lustration). 

The dimensionless surface mass density produced by an SIE 
halo at the location (9[ , 9' 2 ) of a background galaxy is given by 



K(0l,0 2 |0E,e,0o) = — , 



4.2 Simulated fields 

We place a halo, with fiducial parameter values 0e = 10", e = 
0.15, 4>o = and redshift z d = 0.6, at the center of a 1.5' x 
1.5' field. The Einstein radius 0e for a singular isothermal spherical 
halo with velocity dispersion a v is given by 



Dds 
D s ' 



(18) 



Thus for the case of 0e = 10" and z s = 1.45 (implying Dds/ D s = 
0.5), the velocity dispersion is about 840 km/s, corresponding to a 
large group or a cluster. 

We then place 80 background galaxies in the field at ran- 
dom positions. A redshift is assigned to each galaxy according to a 
Gamma distribution with zq = 1/3, 



p(z) = 



2z 3 



exp 



20 



(19) 



(13) 



which peaks at z — 2/3 and has a mean redshift of (z) = 3zo = 1. 
The source density corresponds to space-based observing condi- 
tions. It is compatible with the density of the weak lensing source 
galaxies observed in the COSMOS field (Schrabback et al. 2010). 

The lensing shear and flexion components are computed for 
each galaxy according to the formula given in the previous section 
and using the cosmological weight function defined in Eq. (4). 

We add noise to the shear and flexion signals. The shear noise 
is generated from a Gaussian distribution with a e = 0.2 per shear 
component, and the flexion noise is generated from a Gaussian dis- 
tribution with crjr = 0.03/" and ag = 0.04/" per flexion com- 
ponent. Galaxies with shear absolute value larger than 0.9 are dis- 
carded from the analysis. This corresponds to about a few percent 
of the sample. We further discard, from the flexion analysis: 1) 
galaxies with flexion absolute value larger than 0.5/", since cor- 
rect flexion estimates cannot be obtained in very strongly distorted 
images (see Schneider & Er 2008). 2) galaxies with small flexions 
(|.F| < 0.005/"), which are exceedingly difficult to measure due to 
their low signal-to-noise ratios. 3) galaxies located at high redshift 
(z > 1.0), which are too faint or too small to allow for a reliable 
flexion measurement (Okura et al. 2008). 4) galaxies located at low 
redshift (z < 0.61) are also discarded since they are either below 
the lens redshift or too close to the lens to be efficiently lensed. In 
the end, we discard about 70% of the flexion data in our analysis. 
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Figure 2. Marginalized SIE 68%, 95% and 99% credible regions using shear data only (dashed) and combined shear + flexion data (solid). Left panel: (0E-e) 
plane. Right panel: (<f>o-t) plane. The cross shows the fiducial input model. 



4.3 SIE results 

We perform a likelihood analysis to assess how well the fiducial 
model can be recovered given the noisy shear and the two noisy 
flexion components in the simulated data. 

The theoretical predictions are calculated using Eqs. (15)-(17) 
and the three parameters (9e, £ and <j>o are varied. We use the fol- 
lowing flat priors on the 3 free parameters. Firstly, the orientation 
4>o could in principle be constrained independently by the shape 
of the luminous host object. However, there might be a systematic 
difference between the host orientation and the dark matter halo 
orientation. This misalignment appears to be small for elliptical 
galaxies (Keeton et al. 1998). We assume the polar angle of the ori- 
entation of the halo is distributed in a range smaller than 7r/4 with 
respect to the known orientation of the host. Secondly, we restrict 
the ellipticity to the range between 0.0 and 0.4. Ellipticity can be 
constrained using the morphologies of the lensing hosts. Finally, 
cluster studies using complementary approaches from kinematics 
and X-ray observations should allow us to independently constrain 
the halo mass, or 8e- We allow about 30% uncertainty around the 
input value and restrict (9e to the range [7", 13"] in our analysis. 

Figure 2 shows the resulting credible contours for the SIE halo 
parameters. Both shear and flexion amplitudes increase with halo 
mass and also increase, on average over source locations, with halo 
ellipticity, as shown in Fig. 3. Accordingly, the likelihood analy- 
sis produces an anti-correlated contour, as shown in Fig. 2 (left 
panel). Flexion is more sensitive to a change in the halo elliptic- 
ity than shear, producing tighter constraints. For the SIE model the 
introduction of flexion data produces a gain of a factor of 2 on the 
marginalized errors of ellipticity and mass (see Tab.l). 

There is also a large improvement on the uncertainty of the 
halo orientation when using flexion data, as seen in the right panel 
of Fig. 2. To check the robustness of the result against sample vari- 
ance we made various realizations of the noise and location of the 
background galaxies, having obtained similar results. We also made 
tests using a higher number of background galaxies, obtaining cor- 
respondingly tighter constraints. 

These results assume perfectly known redshifts for the source 
galaxies. We introduce now a redshift uncertainty of 2% on average 
to simulate photo-z errors, which increase with redshift (Bolzonella 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

e 

Figure 3. Amplitude of shear (dashed) and flexion (solid) at a typical loca- 
tion (6[ , 6»2)=(10", 30") as function of the SIE halo ellipticity, normalized 
by the amplitude in the singular isothermal spherical model. 

et al. 2000; Hildebrandt et al. 2010). The marginalized contours for 
SIE parameters are shown in Fig. 4. The impact of the redshift un- 
certainty is significant, due to the strong degeneracy between #e 
and redshift, increasing the la error on the #e estimate by roughly 
a factor of 2. Furthermore, the redshift uncertainty produces a selec- 
tion bias in realizations with excess of low-redshift sources, intro- 
ducing a bias in the #e estimate. Due to the anti-correlation found 
between #e and e in the SIE model, such selection bias will also 
affect the ellipticity estimate. On the other hand, we found no sig- 
nificant effect on the orientation. 



4.4 eNFW model 

In order to test the universality of the result we perform our anal- 
ysis on a more realistic model of halo density profile, adopting 
an eNFW model for the foreground lens halo. The NFW profile 
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SIE 


eNFW 




Shear 


Shear + Flexion 


Shear 


Shear + Flexion 




0.115 ±0.061 


0.164 ±0.035 


0.242 ±0.101 


0.174 ±0.060 


M") 


10.38 ±0.78 


9.75 ± 0.46 












0.256 ± 0.030 


0.242 ± 0.022 


<5M 


15% 


9% 


32% 


26% 



Table 1. lcr errors estimated for the parameters of both models using shear and shear+flexion data. For comparison between SIE and eNFW the uncertainties 
on lens strength parameters were converted into mass uncertainties in the last row. The mass of SIE model stands for the mass within the Einstein radius (#e)> 
and the mass of eNFW model stands for M200. 




0.0 0.1 



0.2 
e 



0.3 



0.4 



Figure 4. SIE 68%, 95% and 99% credible regions using shear and flexion 
data, marginalized over redshift uncertainty. The result of Fig. 2 is shown 
for comparison (solid contours). The cross shows the fiducial input model. 



(Navarro et al. 1996, 1997) is widely used to model the halo of 
galaxy clusters. The dimensionless surface mass density of a spher- 
ical NFW halo is written as (Bartelmann 1996; Bacon et al. 2006) 

/(*) 



k(x) — 2n s 



r 



(20) 



where x is the dimensionless radius, the radius r normalized by the 
scaling radius r B (x = r/r s ), and the function f(x) is given by 



fix) = < 



zarctan 



x- 1 
x + 1 



(x > 1); 



(z = l); 



(21) 



zarctanh 



\+x 



(X<1). 



The physical properties of the halo are contained in the parameter 
k s = pcrit7* s A c /Ecr, where A c is the dimensionless characteristic 
density. The halo mass is defined as M200 = 4/3p C rit7iT2nn, where 
r200 — r B c and c is the concentration parameter (see the appendix 
in Navarro et al. 1997). 

For an elliptical halo the dimensionless radius at a point 
(61 ,62), defined along the axes of the halo, becomes 

p D d 



1 — e r s 

The halo ellipticity is defined by e = 



(0a 



(22) 

0b)/(0 a + 0b) (or 



equivalently by the axial ratio / = 6 b /9 a ), where 8 a , 9 b are the 
major and minor axes. 

The lensing properties of the eNFW halo can be calculated nu- 
merically given an arbitrary normalized halo convergence (Keeton 
2001; Hawken & Bridle 2009) as follows. 

The shear and flexion are second and third order derivatives of 
the lensing potential rp: 



7 

g 



(ipll - V>22) + iV"i2; 

[ipm + V»i22 + i(V>ii2 + ^222)] ; 

[iprii - 3^122 + i(3V>ii2 - ^222)] , 



(23) 
(24) 
(25) 



where subscripts denote partial differentiation. The second and 
third order derivatives of the lensing potential, for an eNFW dark 
matter halo at a position (x, y) on the image plane, are given by 



V11 = 2fx 2 K + fJo; 

V22 = 2fy 2 K + fJ 1 ; 

V>i2 = 2fxyKr, 

tpiu = 6fxK + 4fx 3 L ; 

1P222 = 6fyK 2 +4fy 3 L 3 ; 

V>ii2 = 2fyK 1 +4:fx 2 yL 1 ; 

V122 = 2fxK 1 +4fy 2 xL 2 . 



Here 

Jn(x,y) 

K n (x,y) 
L n (x,y) 



- f 

Jo 

- f 

Jo 



K(C(u) 2 )du 



[1 - (1 - />r +i/a ' 

lift' (C( w ) 2 ) du 

[1 - (1 - / 2 )u]"+ 1 /2 ' 

r 1 u 2 k" (ij(u) 2 )du 

Jo [1 - (1 - rw+w 



(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 



(33) 
(34) 
(35) 



are one-dimensional integrals, where / is the axis ratio of the lens 
and k and k" are the first and second order derivatives of the con- 
vergence, e.g. k'(£ 2 ) = dK(£ 2 )/d(£ 2 ). The convergence k is writ- 
ten as a function of the ellipse coordinate £ (u) given by 



= u [x + 



y 



1-(1-P)u 



(36) 



Notice here we need to use x 2 instead of x as the independent vari- 
able in k(x). 
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4.5 eNFW results 

In the eNFW model, we consider 3 free parameters: (k s , e, (j>). For 
our fiducial halo we use a halo mass of M200 = 1-8 x 1O 14 M , a 
concentration parameter of c = 7.2, and place the halo at redshift 
z d = 0.6. This implies k s — 0.227. For the ellipticity we choose 
e = 0.15 and for position angle <po = 0. 

The theoretical predictions are calculated numerically and the 
three parameters « s , e and 4>o are varied. Similar to the SIE anal- 
ysis, we place 80 background galaxies at random positions in the 
same field, and randomize their shear and flexion values. The same 
filter to the data is employed to discard the unmeasurable data in 
the analysis. We perform a likelihood analysis restricting cj>o to the 
range [— 7r/4, it/4], using the range [0.15, 0.3] for k 3 and [0, 0.4] 
for the ellipticity e. 

Figure 5 shows the resulting credible intervals for the eNFW 
halo parameters. Overall, the constraints are looser than in the SIE 
model, which has a higher signal-to-noise ratio. The effective num- 
ber of flexed background images available for the analysis after 
discarding is lower than that in the SIE analysis by about 15%. 

The constraints from shear information alone, using priors 
similar to the ones used in the SIE analysis, are dominated by the 
priors. Hence, shear data do not add much information to observa- 
tions of the luminous host and complementary probes of mass (e.g. 
X-ray). 

Ellipticity and lens strength are now positively correlated in 
the shear signal due to the decreasing shear amplitude with ellip- 
ticity, as shown in Fig. 6. The addition of flexion data decreases the 
correlation. Once again, flexion is more sensitive to the change in 
the halo ellipticity than shear. Indeed, the tangential shear mainly 
depends on the lens strength (9e or n s ) and depends weakly on 
the halo ellipticity, while the cross shear component is independent 
of ellipticity. On the other hand, both flexion components depend 
on the halo ellipticity and orientation <po. This allows for a tighter 
constraint when including flexion information. In Table 1 the shear 
and flexion marginalized constraints on the ellipticity and strength 
parameters are roughly a factor of 1.5 tighter than the correspond- 
ing shear constraints. Notice that the effective gain of using flexion 
is larger than this factor. Here the shear constraint, contrary to the 



V 


1 1 1 
flexion 




shear 


- eNFW halo 
1 1 1 


1 1 1 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



e 

Figure 6. Amplitude of shear (dashed) and flexion (solid) at a typical lo- 
cation (e[,e' 2 ) = (10", 30") as a function of the eNFW halo ellipticity, 
normalized by the amplitude in the spherical NFW model. 

combined one, is dominated by the prior; with a wider prior range 
the shear constraint would be looser. 

We also stress that the eNFW model constraints are not 
marginalized on all halo parameters, since the scaling radius r s is 
kept fixed in the analysis. Finally, we checked the results against 
sample variance making various realizations of noise and location 
of the background galaxies, similar results are found. 



5 CONCLUSIONS 

In this paper, we study the potential of weak lensing flexion in the 
study of galaxy cluster haloes. We use mock data including shear, 
T flexion, Q flexion and redshift information. We find that the in- 
clusion of flexion significantly improves the estimate of foreground 
haloes parameters, although the details are model-dependent. In 
particular, in the case of a SIE halo, the presence of a mass- 
ellipticity anti-correlation implies that analyses where the halo is 
incorrectly assumed to be spherical will overestimate the halo mass. 
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On the other hand, we did not find significant correlation between 
the halo mass and ellipticity in the eNFW model. 

The noise in the mock data is determined by the dispersion 
of the intrinsic shear and flexion distributions, and by the density 
of background galaxies. After applying stringent cuts in the data, 
we are left with a galaxy density of roughly 10 arcmin~ 2 . Our ap- 
proach assumes that the flexion estimators are linear in the flexion 
observables and the point spread function can be removed without 
producing a bias. In reality, the noise of flexion estimators can be 
complicated and non-Gaussian. An accurate study of flexion noise 
is important to evaluate the estimated error. 

The analysis considers a single cluster halo. This approach is 
not possible if the number density of background galaxies is low. In 
that case, stacking of several halo fields can be used to increase the 
number of background images and thus the signal-to-noise. That 
approach requires the alignment of the major axis of various fore- 
ground galaxies and selecting haloes with similar properties, for ex- 
ample similar shapes of their central galaxies. Such stacking analy- 
sis can constrain the halo shapes more tightly, and as a function of 
other halo properties, e.g. mass. 

Our results emphasize that a combined weak lensing analysis 
will be a useful technique for precise measurements of the proper- 
ties of galaxy or cluster haloes from future weak lensing surveys, 
such as EUCLID. 
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